Model independent constraints on contact interactions 
from polarized Bhabha scattering at LC. 



A. Pankov a ' 6 and N. Paver a 

a Dipartimento di Fisica Teorica, Universita di Trieste and 
Istituto Nazionale di Fisica Nucleare, Sezione di Trieste, Trieste, Italy 
b The Abdus Salam International Centre for Theoretical Physics, Trieste, Italy 



Abstract 

We discuss electron-electron contact-interaction searches in the process e + e~ — > e + e~ at 
a future e + e~ Linear Collider with cm. energy y/s = 0.5 TeV and with both beams 
longitudinally polarized. Our analysis is based on the measurement of the polarized dif- 
ferential cross sections, and allows to simultaneously take into account the general set of 
contact interaction couplings as independent, non-zero, parameters thus avoiding simpli- 
fying, model-dependent, assumptions. We evaluate the corresponding model-independent 
constraints on the coupling constants, emphasizing the role of beam polarization, and make 
a comparison with the case of e + e~ — > fi + . 



1 Introduction 



Contact interaction Lagrangians (CI) provide an effective framework to account for the 
phenomenological effects of new dynamics characterized by extremely high intrinsic mass 
scales A, at the 'low' energies ^/s <C A attainable at current particle accelerators. Typically, 
'low energy' manifestations of quark and lepton substructure would occur via four-fermion 
quark and lepton contact interactions, induced by exchanges of very heavy sub-constituent 
bound states with mass of the order of A pQ . However, this concept is quite general, and 
is adequate to the description of a wide class of quark and lepton cross sections governed 
by exchanges of objects associated to new gauge interactions such as, e.g., heavy neutral 
gauge bosons and leptoquarks, in a regime where such masses are much larger than the 
Mandelstam variables of the process. 

The explicit parameterization of the four-fermion quark and lepton contact interactions 
is, a priori, somewhat arbitrary. In general, it must respect SU(3) x SU(2) x £7(1) symme- 
try, because the new dynamics are active well-beyond the electroweak scale. Furthermore, 
usually one limits to the lowest dimensional operators, D = 6 being the minimum, and 
neglects higher dimensional operators that are suppressed by higher powers of 1/A 2 and 
therefore are expected to give negligible effects. 

In this note, we consider the effects of the flavor- diagonal, helicity conserving, eeff 
contact-interaction effective Lagrangian [2j 

£ci = : — ; Y] 9es (ea^) (frffj) , (1) 
l + Oef r: 

J 1,3 

in the Bhabha scattering process 

e + + e~ -y e + + e" (2) 

at an e + e~ Linear Collider (LC) with cm. energy y/s = 0.5 TeV and polarized electron 
and positron beams. In Eq. (0): i, j — L, R denote left- or right-handed fermion helicities, 
/ indicates the fermion species, so that 8 e f = 1 for the process (j2j) under consideration, 
and the CI coupling constants are parameterized in terms of corresponding mass scales as 
6ij = rjij/Afj. Actually, one assumes g1 s = 4tt to account for the fact that the interaction 
would become strong at y/s ~ A and, by convention, |^-| = ±1 or = 0. This leaves the 
energy scales Ajj as free, a priori independent parameters. 

Clearly, at s A?-, the Lagrangian ((TJ) can only contribute virtual effects, to be sought- 
for as very small deviations of the measured observables from the Standard Model (SM) 
predictions. The relative size of such effects is expected to be of order s/aA 2 , with a 
the SM coupling (essentially, the fine structure constant) and, therefore, very high collider 
energies and luminosities are required for this kind of searches. In practice, the constraints 
and the attainable reach on the CI couplings can be numerically assessed by comparing the 
theoretical deviations with the foreseen experimental uncertainties on the cross sections. 

For the case of the Bhabha process (J2J), the effective Lagrangian Cci in Eq- O envisages 
the existence of six individual, and independent, CI models, contributing to individual he- 
licity amplitudes or combinations of them with a priori free, and nonvanishing, coefficients 
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(basically, £rr and €lr combined with the ± signs). Correspondingly, in principle the 
most general, and model- independent, analysis of the data must account for the situation 
where all four-fermion effective couplings defined in Eq. (£TJ are simultaneously allowed in 
the expression for the cross section. Potentially, the different CI couplings may interfere 
and substantially weaken the bounds. Indeed, although the different helicity amplitudes 
by themselves do not interfere, the deviations from the SM could be positive for one helic- 
ity amplitude and negative for another, so that accidental cancellation might occur in the 
sought-for deviations from the SM predictions for the relevant observables. 

The simplest attitude is to assume non-zero values for only one of the couplings (or one 
specific combination of them) at a time, with all others zero, and this leads to tests of the 
specific models mentioned above. Also, in many cases, global analyses combining data from 
different experiments relevant to the considered type of coupling are performed. Current 
lower bounds on the corresponding A's obtained along this line from recent analyses of 
e + e~ — * // at LEP, that include Bhabha scattering, are in the range 8-20 TeV and are 
found to substantially depend on the considered one-parameter scenario [31 H] . Examples 
of results for the eeff couplings for the different fermion species in from analyses of 
different kinds of processes and experiments, can be found, e.g., in Refs. [SHU]- 

It should be highly desirable to apply a more general (and model-independent) ap- 
proach to the analysis of the experimental data, that allows to simultaneously include all 
terms of Eq. as independent, non vanishing free parameters, and yet to derive separate 
constraints (or exclusion regions) on the values of the CI coupling constants, free from 
potential weakening due to accidental cancellations. 

Such an analysis is feasible with initial beam longitudinal polarization, a possibility 
envisaged at the LC ^Hj- This allows to extract the individual helicity cross sections 
from suitable combinations of measurable polarized cross sections and, consequently, to 
disentangle the constraints on the corresponding CI constants e^, see, e.g., Refs. [TTHrH] . 
In what follows, we wish to complement the model-independent analysis of e + e~ — > // 
with f e,t given in Refs. fTjm], with a discussion of the role of the polarized differential 
cross sections measurable at the LC in the derivation of model-independent bounds on the 
three independent four-electron contact interactions relevant to the Bhabha process 00 

Specifically, in Sect. 2 we introduce the observables being considered, and in Sect. 3 we 
present a numerical analysis and the assessment of the attainable reach on the CI couplings, 
resulting from a \ 2 procedure that accounts for the expected experimental uncertainties 
(statistical and systematical ones). Finally, Sect. 4 contains a comparison with the results 
previously obtained for e + e~ — > [i + and some conclusive remarks. 

2 Polarized observables 

With P~ and P + the longitudinal polarization of the electron and positron beams, respec- 
tively, and 9 the angle between the incoming and the outgoing electrons in the cm. frame, 

Notice that, in general, for e + e~ — * // with / ^ e there are four independent CI couplings, apart 
from the ± possibility, so that in the present case of process there is one free parameter less. 
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the differential cross section of process (J2J) at lowest order, including 7 and Z exchanges 
both in the s and t channels and the contact interaction ((TJ), can be written in the following 
form [lU [121 EH]: 

MP',P + ) p - p+) (1 p - p+) (p+ _ p-x (3) 

In Eq. ©: 



do"i 


7ra 2 


dcos^ 


4s 


dcx 2 


71a 2 


dcos^ 


4s 


dap 


na 2 


dcos# 


4s 



COS I 



[A + (l + cos#) 2 + A_( 
Aj(l + cos0) 2 , (4) 



with 



|2 

1 



^o(s,t) = (7) |l + #r#lXz(*) + -clr| 
A+M = 1 1 1 + I + (xz{s) + S - X z{t)) + 2 ^e LL | 2 

5 2 

A_(M) = 1 + fl'Rfl , LXz(-5) + -e LR , 

^JM) = 1 1 1 + I + ^£ (x^(s) + |xz(*)) + 2 ^ e L L | 2 

- + J + ^r (xz(s) + yXz(t)) +2^e RR | 2 . (5) 

Here: a is the fine structure constant; t = — s(l— cos#)/2 and Xz(s) = s/(s— M^+iMzT z), 
Xz(t) = t/(t — M|) represent the Z propagator in the s and t channels, respectively, with 
M z and the mass and width of the Z; g R = tan#^, 9h — ~ cot 2 6 W are the SM right- 
and left-handed electron couplings of the Z, with 8\y the electroweak mixing angle. 

With both beams polarized, the polarization of each beam can be changed on a pulse 
by pulse basis. This would allow the separate measurement of the polarized cross sections 

for each of the four polarization configurations ++, , H — and — h, corresponding to the 

four sets of beam polarizations (P~, P + ) = (Pi, P2), (—Pi, —P2), (Pi, —P2) and (—-Pi, P2), 
respectively, with < Pi 2 < 1. Specifically, with the simplifying notation da = da/dcos6: 

da ++ = da(P 1 ,P 2 ) = (l-P 1 P2)dai + (l + PiP 2 )da2 + (P2-Pi)dap, 

da__ = da(-Pi,-P 2 ) = (l-PiP2)da 1 + (l + P 1 P 2 )da 2 -(P2-P 1 )dap, 

da+_ = da(Pi,-P 2 ) = (l + P 1 P 2 )dai + (l-PiP 2 )da 2 -(P 2 + P 1 )da P , 

da_ + = da(-Pi,P 2 ) = (l + PiP 2 )da 1 + (l-P 1 P 2 )da 2 + (P 2 + P 1 )d(Tp. (6) 
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To extract from the measured polarized cross sections the values of d<7i, d<T2 and do>, that 
carry the information on the CI couplings, one has to invert the system of equations 
The solution reads: 

+ da_+) , 
+ da_+) , 
" da_) • (7) 

Notice that the equations in (jUJ) are not all linearly independent, and that not only Pi ^ 

and P2 7^ 0, but also P\ 7^ P 2 is needed to obtain dap via da ++ and da As one can 

see from Eqs.(JH) and a 2 depends on only one contact interaction parameter (€lr), 
ap is two-parameter dependent (€rr and and a\ depends on all three parameters. 
Therefore, the derivation of the model-independent constraints on the CI couplings requires 
the combination of all polarized observables of Eq. ©• In this regard, to emphasize the role 
of polarization, one can observe from Eqs. ©-© that, in the unpolarized case Pi = P 2 = 
where only a\ and <7 2 appear, the interference of the €lr term with the SM amplitude in Aq 
and A_ has opposite signs, leading to a partial cancellation for —t ~ s. Consequently, as 
anticipated in Sect. 1, one expects the unpolarized cross section to have reduced sensitivity 
to e L R. Conversely, £lr is directly accessible from d<r 2 , via polarized cross sections as in 
Eq. (JJJ). Also, considering that numerically g\ = g^, the parameters €ll and e^R contribute 
to the unpolarized cross section through A + with equal coefficients, so that, in general, 
only correlations of the form + £rr| < const, and not finite allowed regions, could be 
derived in this case. 

To make contact to the experiment we take P\ = 0.8 and P2 = 0.6, and impose a cut 
in the forward and backward directions. Specifically, we consider the cut angular range 
I cos 6*| < 0.9 and divide it into nine equal-size bins of width Az = 0.2 (z = cos 9). We also 
introduce the experimental efficiency, e, for detecting the final e + e~ pair and, according to 
the LEP2 experience, e = 0.9 is assumed. 

We then define the four, directly measurable, event rates integrated over each bin: 

N ++ , AL„, iV + _, AL + , (8) 

and (a/3 = ++, etc.): 

= -Ante / (da aP /dz)dz. (9) 

4 J bin 

In Eq. (JDJ), £j nt is the time-integrated luminosity, which is assumed to be equally di- 
vided among the four combinations of electron and positron beams polarization defined in 
Eqs. ®. 

In Fig. 1, the bin- integrated angular distributions of iV'j"" and iV^ 1 " in the SM at ^/s = 
500 GeV and C mt = 50 fbT 1 are presented as histograms. Here, the SM cross sections have 

4 
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2(Pi + P 2 ) 



(da, 



da. 



2(p 2 -n: 



(da, 
(da, 
(da, 



been evaluated by means of the effective Born approximation [TTH |2"U] . The typical forward 
peak, dominated by the t-channel photon pole, dramatically shows up, and determines a 
really large statistics available in the region of small t. The cos 9 distributions for the other 
polarization configurations in Q are similar and, therefore, we do not represent them here. 
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Figure 1: Bin-integrated angular distributions of (solid line) and (dashed line), 
Eq.@, in the SMat^ = 500 GeV and £ int = 50 fb -1 . 



The next step is to define the relative deviations of the cross sections a%, a 2 and op 
from the SM predictions, due to the contact interaction. In general, for such deviations, 
we use the notation: 



A 



o 



Q{SM + CI) - Q{SM) 
0{SM) ' 



(10) 



To get an illustration of the effect of the contact interactions on the observables (JZJ) under 
consideration, we show in Fig. 2a and Fig. 2b the angular distributions of the relative 
deviations of d<7i and da 2 , taking as examples the values of £ int and Ay indicated in the 
captions. The SM predictions are evaluated in the same, effective Born, approximation 
as in Fig. 1. The deviations are then compared to the expected statistical uncertainties, 
represented by the vertical bars. Fig. 2a shows that dcxi is sensitive to contact interactions 
in the forward region, where the ratio of the 'signal' to the statistical uncertainty increases. 
Also, it indicates that, for the chosen values of the cm. energy ^/s and £ int , the reach on 
will be substantially larger than 30 TeV. Conversely, Fig. 2b shows that the sensitivity of 
d<72 is almost independent on the chosen kinematical range in cos#, leading to a really high 
sensitivity of this observable to €lr, and to corresponding lower bounds on Alr potentially 
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larger than 50 TeV. We now proceed to the analysis of the bounds on the contact interaction 
couplings. 
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Figure 2: The angular distributions for the relative deviations A ai (a) at A RR =30 TeV 
(solid line) and 50 TeV (dashed line), and for A CT2 (b) at A LR =40 TeV (dashed line), 
50 TeV (solid line), 70 TeV (dot-dashed line). The curves above (below) the horizontal 
line correspond to negative (positive) interference between contact interaction and SM 
amplitude. The error bars show the expected statistical error at £ int = 50 fb~ . 



3 Numerical analysis and constraints on CI couplings 

To assess the sensitivity to the compositeness scale we assume the data to be well described 
by the SM predictions (e^ = 0), i.e., that no deviation is observed within the foreseen 
experimental accuracy, and perform a x 2 analysis of the cos# angular distribution. For 
each of the observable cross sections, the x 2 distribution is defined as the sum over the 
above mentioned nine equal-size cos# bins: 

bins ^ ' 

where O = 01, a 2 , o> and a hm = J Un (da/dz)dz. In Eq. (ITTj) . AO represents the deviation 
from the SM prediction, AO = 0(SM + CI) — O(SM), which can be easily expressed 
in terms of the CI couplings by using Eqs. (jHJ), and SO is the expected experimental 
uncertainty, that combines the statistical and the systematic ones. 

In the following analysis, the theoretical expectations for the polarized cross sections 
are evaluated by using the program TOPAZ0 [2B 122] , adapted to the present discussion, 
with m top = 175 GeV and M# = 120 GeV. For electron-positron final states, a cut on 
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the acollinearity angle between electron and positron, # aco i < 10°, is applied to select 
non-radiative events. 

Concerning the numerical inputs and assumptions used in the estimate of SO, to as- 
sess the role of statistics we vary £ int from 50 to 500 fb" 1 (a quarter of the total run- 
ning time for each polarization configuration). As for the systematic uncertainty, we take 
5C int /C int = 0.5%, 5e/e = 0.5% and, regarding the electron and positron degrees of polar- 
ization, 5P 1 /P 1 = 5P 2 /P 2 = 0.5 %. 

For example, in the case of a 2 , the deviation from the SM prediction, largely dominated 
for y/s << Alr by the interference term, can be represented as: 

Aa 2 bin = af^SM + CI) - af*(SM) ~ 2nase LR [ (1 + g^ X z{t)) , (12) 

J bin t 

and, combining uncertainties in quadrature, the uncertainty on a 2 , indirectly determined 
via the measured a ++ , a , and as in Eq. ((7|), can be expressed as 



(Sa 2 f 
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■(13) 



with 



5a, 



a/3 



a, 



a/3 



N, 



a/3 



5C 



hit 



c 



mi 



(14) 



Analogous expressions hold for the uncertainties relevant to a\ and ap. 

As a criterion to constrain the allowed values of the contact interaction parameters 
by the non-observation of the corresponding deviations, we impose x 2 < Xcl' where the 
actual value of Xcl specifies the desired 'confidence' level. We take the values Xcl =3-84 
and 5.99 for 95% C.L. for a one- and a two-parameter fit, respectively. 

We begin with the presentation of the numerical results for €lr. In this case, the relevant 
cross section a 2 depends on clr only, see Eqs. (HJ) and (J3J) and, therefore, the constraints 
on that parameter are determined from a one-parameter fit. The model-independent, 
discovery reach expected at the LC for the corresponding mass scale A L r is represented, 
as a function of the integrated luminosity Ant, by the solid line in Fig. 3. As expected, 
the highest luminosity determines the strongest constraints on the CI couplingso Fig. 3 
dramatically shows the really high sensitivity of a 2 , such that the discovery limits on 
A L r are the highest, compared to the A RR and A L l cases, and can be as large as 110 up 
to 170 times the total cm. energy. Actually, the Alr limits are large enough that the 
approximation used in Eq. (fT2|) is good for all values of t, so that the derived limits are 
independent of the sign of the 6lr. 

2 Such increase with luminosity is somewhat slower than expected from the scaling law A ~ (sCint) 1 ^ 4 , 
as the effect of the systematic uncertainties competing with the statistical ones. 
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Figure 3: Reach in Ay at 95% CL. vs. integrated luminosity C mt obtained from the 
model-independent analysis for e + + e~ — > e + + e~ at y/s = 0.5 TeV, \P~\ = 0.8 and 
|P + | = 0.6, A L r (solid line), A L l (dashed line), A RR (dot-dash line). 



Since op simultaneously depends on the two independent CI couplings £rr and 6ll, 
a two-parameter analysis is needed in this case. Here, the terms quadratic in eix and 
£rr largely cancel leaving the remaining interference to dominate the relevant deviations 
from the SM, and consequently the resulting constraint has the form of a straight band, 
as depicted in Fig. 4a. Indeed, such a band represents a correlation between the two 
parameters, rather than a bound around the SM value €ll = crr = 0. 

In order to get a restricted allowed region around zero, one can combine the above 
mentioned band with the exclusion region obtained from a±. However, since the latter 
depends on all three contact interaction parameters, see Eqs. (@J and ©, to set constraints 
in the (€rr, €ll) plane requires the combination of the o"i-bounds with the limits on €lr 
derived above from <r 2 . The bound in the (crr, £ll) resulting from this procedure is shown 
in Fig. 4a and, finally, the shaded ellipse determined by the combination with the band 
determined by op represents the restricted allowed region around the SM point e^- = 0. 
With reference to Eq. (fTTj) . for the x 2 analysis this amounts to the consideration of the 
combined x 2 (°"i) + X 2 ( cj p)- Fig- 4b is essentially a magnification of the shaded region of 
Fig. 4a, and represents the model-independent limits on e LL and £rr attainable at the 
considered LC, for two possible values of the integrated luminosity. These bounds are 
translated into the model-independent reach on the mass scale parameters All and Arr, 
represented as a function of luminosity in Fig. 3. The fact that such bounds are substan- 
tially lower than those for Alr reflects that a combined two-parameter \ 2 analysis must be 
used. In this regard, the calculation presented here indicates that not only polarization, 
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but also combinations of measurements of polarized observables are necessary to obtain 
model-independent bounds on the CI couplings. 




Figure 4: (a) Allowed areas at 95% CL. on electron contact interaction parameters in 
the planes (crr^ll), obtained from o\ and <jp at y/s = 500 GeV and C- mt = 50 fb _1 (a); 
(b) Combined allowed regions at 95% C.L. obtained from a\ and op at y/s = 500 GeV, 
Ant = 50 fb~ (outer ellipse) and 500 fb _1 (inner ellipse). 



The crosses in Fig. 4b represent the model-dependent constraints obtainable by taking 
only one non-zero parameter at a time, either eix or £rr, instead of the two simultaneously 
non-zero and independent as in the analysis discussed above. Similar to the inner and 
outer ellipses, the shorter and longer arms of the crosses refer to integrated luminosity 
Ant — 50 fb _1 and 500 fb _1 , respectively. One can note from Fig. 4b that the 'single- 
parameter' constraints on the individual CI parameters £rr and ejx are numerically more 
stringent, as compared to the model-independent ones. Essentially, this is a reflection of 
the smaller value of the critical x 2 , Xcl = 3-84 corresponding to 95% C.L. with a one- 
parameter fit. 



4 Concluding remarks 

In the previous sections we have derived limits on the contact interactions relevant to 
Bhabha scattering by a model-independent analysis that allows to simultaneously account 
for all independent couplings as non-vanishing free parameters. The results for the lower 
bounds on the corresponding mass scales A range, depending on the luminosity, from 
essentially 38 to 50 TeV for the LL and RR cases, and from 54 to 84 TeV for the LR case. 
The comparison with the numerical results relevant to the e + e~ — > channel, derived 

from a similar analysis [18 , is shown in Table 1. 

The table shows that for All and Arr the restrictions from e + e~ — > and e + e~ — > 

e + e~ are qualitatively comparable. Instead, the sensitivity to Alr, and the corresponding 
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Table 1: Reach in A^- at 95% C.L., from the model-independent analysis performed for 
e+e" -> [i+fi- and e + e~, at E c , m , = 0.5 TeV, C mt = 50 ftT 1 and 500 fb" 1 , \P~\ = 0.8 and 
IP+I = 0.6. 



process 


Ant 


A LL 


Arr 


A L r 


Arl 
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TeV 


TeV 


TeV 


TeV 




50 


35 


35 


31 


31 


e + e" — > 


500 


47 


49 


51 


52 




50 


38 


36 


54 




e + e~ — > e + e" 


500 


51 


49 


84 





lower bound, is dramatically higher in the case of Bhabha scattering. In this regard, this 
is the consequence of the initial beams longitudinal polarization that allows, by measuring 
suitable combinations of polarized cross sections, to directly disentangle the coupling €lr. 
Indeed, as previously observed, in general without polarization only correlations among 
couplings, rather that finite allowed regions, could be derived or, alternatively, a one- 
parameter analysis testing individual models can be performed. 

As an example of application of the obtained results to a possible source of contact 
interactions, we may consider the sneutrino parameters (mass and Yukawa coupling 
A) envisaged by supersymmetric theories with 7£-parity breaking. In this case, sneutrino 
exchange affects only those helicity amplitudes with non-diagonal chiral indices, so that 
Alr is the relevant mass scale [2S1I2I!- Qualitatively, without entering into a detailed and 
more complex analysis, one can expect typical bounds on m^/A ~ A L r/a/87t ~ 11 to 17 
TeV corresponding to A LR « 54 TeV and 84 TeV (Fig. 3) at C mt = 50 fb" 1 and 500 fb" 1 , 
respectively. 
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